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Abstract. 

The concept of thermal ratchets is extended to the system governed by quan- 
tum mechanics. We study a tight-binding model with an asymmetric periodic 
potential contacting with a heat bath under an external oscillating field as a 
specific example of quantum ratchet. Dynamics of a density operator of this 
system is studied numerically by using the quantum Liouville equation. Finite 
net current is found in the non-equilibrium steady state. The direction of the 
current varies with parameters, in contrast with the classical thermal ratchets. 
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Brownian particles under spatially asymmetric potential attract much at- 
tention recently, which has been shown to produce nonzero net current un- 
der certain conditions^ |2|. This type of the system is called thermal ratchet 
because particles are subject to the thermal noise. In the context of biol- 
ogy, thermal ratchets has been discussed as a possible mechanism of biological 
motors [|| U |[ |). They also can be a mechanism of actuators for micro- or 
nano-machines. In the field of physics, it is studied in connection with the noise 
induced transportation § 

For nonzero net current to appear in thermal ratchets both the spatial and 
time reversal symmetries need to be broken. Otherwise, no finite net current 
is expected, because of the Curie's principle fl6||, which states that asymmetric 
phenomena is caused by asymmetry (when no symmetry breaks spontaneously). 
Hence both (1) existence of an spatially asymmetric potential and (2) a contact 
with a heat bath, which breaks the time reversal symmetry of the system by 
energy dissipation, are required for the realization of thermal ratchets. These 
two requirements alone, however, can not produce finite net current due to the 
second law of thermodynamics. In order to make system nonequilibrium, and 
thus to allow steady current, we require (3) existence of an additional noise or a 
time correlated force. A spatially unbiased noise or force is assumed hereafter, 
because otherwise an occurrence of directional motion of the particle is a trivial 
consequence of the asymmetry. We expect a directional particle motion in a 
non-equilibrium steady state, if and only if above three conditions are fulfilled. 

The models of thermal ratchets discussed so far in the literatures are based 
on classical mechanics. But there are situations in which quantum effects such 
as quantum tunneling play important roles in the particle motion; for instance, 
in small systems like mesoscopic devices at low temperature. Here we propose a 
new concept, quantum ratchets, which are ratchet systems governed by quantum 
mechanics. Quantum ratchets give rise to interesting questions about nonequi- 
librium quantum dynamics in the presence of interaction with environments. 
Moreover, by studying quantum ratchet systems, we would be able to explore 
the new paradigm of quantum thermodynamics such as a quantum engine or a 
quantum energy transducer. In addition, from the technological point of view, 
quantum ratchets have a potential of being used as a new quantum electronic de- 
vices or quantum motors for nano-machines in future. In this letter, we propose 
one possible realization of quantum ratchets and show that finite net current 
appear in that system, which originates from purely quantum mechanical effects. 

The model investigated in this letter is a tight-binding model. Its Hamilto- 
nian is given as follows: 

H to tai{t) = H sys (t)+-f H m t+H B , (1) 

H S y S 

(*) = 12 i\ n )( n + M + H(n-l\ + (V n mod N + F n (t))\n)(n\}, (2) 
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H B = ^hcj a (aia a + ^j, (3) 

■Hint = X! {^^+ll n )( n + !| + 0n-l,n|™)(™ - 1| 
+ 72 K mod AT l«-)H } (Oa + «a) 

= Hz^iaa + at). (4) 

a 

H total (t) is Hamiltonian of the total system consisting of the ratchet part and the 
heat bath. The ratchet part is described by H sys (i) . The localized state at site 
n is represented by \n). V n mo d n represents an asymmetric periodic potential 
with period N and F n (t) an external time dependent field. We measure energy 
in unit of the hopping coefficient. The heat bath, which is described by Hb, 
is a collection of harmonic oscillators, each of which has a frequency uj a and is 
described by the creation(annihilation) operator a^(a a ). The ratchet part and 
the heat bath interact with each other with a coupling constant 70 via Hi nt , 
which is factorized into the following two parts: Hj and X)a( a a + a a)i which 
operate on the Hilbert space of the ratchet part and the heat bath, respectively. 
In Hi the function g n . m is defined as 1 + "fi(V n mo d n + ^mmodiv) with the 
coupling parameter 71; 72 is another coupling parameter. This interaction is a 
general form of real symmetric tridiagonal interactions up to the first order of 
V„ (We omit 0-th order of the diagonal part, since it is a multiple of an identity 
operator and does not contribute to the energy dissipation). We note that the 
Hamiltonian Htotai (i) is Hermitian, because only the diagonal part of Htotai (t) 
is spatially asymmetric. 

Let us start with the quantum Liouville equation for the density operator of 
the total system <j{t), 

^- = ~[H total (t),a(t)}. (5) 

The density operator of the ratchet part p(t) is obtained by p(t) = Tr B<j(t), 
where Tr b means the trace operation with respect to the Hilbert space of the 
heat bathjl7|, Assuming that (1) the total density operator at initial time, 
t = 0, is decomposed into the ratchet part and the heat bath part, i.e., er(0) = 
p(0)pB, (2) the heat bath is at its thermal equilibrium at temperature Tb at h, so 
that p B = exp (-/3H B ) /Tr B exp (~/3H B ) , (3 = l/k B T bath , and (3) the coupling 
constant 70 is sufficiently small, we get the following equation for p(t) after 
standard calculations: 

^ = ~[H sys (t),P(t)} 
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- [H I ,p(t-a)*' t -'H t I ' t -'M-e)), (6) 

where &(t) is the autocorrelation function of the operator £(£) = exp{(i/?i)iJst} 
Y, a (a a + al)e-xp{-(i/h)H B t}, i.e., §(t) = Tr B {p B £(t)£(0)}. We introduced 
the notation A* ,T for an arbitrary operator A as follows: 

A l ' T = *T cxp | - ~ ^ tf sys (f )df | AT* exp j - jf ff sys (f )df | , (7) 

where T, X indicate time ordering product from the right to the left, and one 
from the left to the right, respectively. Equation (||) is exact up to order 7q. 

Since we have assumed the parameter 70 very small, the characteristic time 
of the relaxation of pit) is very long. Moreover, if the heat bath is sufficiently 
large, the function <!>(£) vanishes quickly. In such situation the integrand near 
s ~ mainly contributes to the integral. In addition, assuming the external field 
changes more slowly than the characteristic time of we can approximate 
eq. (§) as 

^ = -~[H sys it),p(t)} - ^([Hj, [K(t),p(t)]]+i[Bi, [H(t),p(t)] + ]), (8) 
where [A, B] + = AB + BA. Here we introduced new operators, 

K(t) = j dsexp(-^H sys (t)s^H I exp(^H sys (t)s^Re^(s) (9) 

Hit) = J dsexp^-^H sys (t)s^H I expl.^H sys (t)s^Im$(s). (10) 

In what follows, we study the dynamics of the system with an asymmetric 
parabolic potential of period 5; i.e., V mod 5 = 0, Vi mo d 5 = V 9 > ^2 mod 5 = 4/9, 
V3 mod 5 = 1, and V4 mod 5 = 1/9. Thus an unit cell consists of five sites. 
For modeling of the heat bath, we use the spectral density of oj which pro- 
duces the standard Ohmic dissipation with the cutoff frequency 1/A[[l9| |2?j, 
then we can replace a summation over modes by an integration, i.e., ^2 a — > 
Jg 00 duiu) exp(— Xui) in the evaluation of $(s). The external field is simply chosen 
as a periodically on-off type ("flushing") field of the period T|| [To|| ; 

p m _ / -V n mo d n for i mod T > T/2 , . 

nlj ~\0 for t mod T < T/2. [ > 

A current is defined through the conservation of the probability; taking the 
diagonal element of the equation of motion (||), we can define the local cur- 
rent J n (t) at the site n, by using the conservation relation for the probability 
{n\(dp{t)/dt)\n) + J n+ i(t) - J n (t) = 0, as follows: 



Jn(t) = -hm(n\p(t)\n - 1) + ^§-g n -i..Jle{n\Wit)\n - 1), (12) 

fl 
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where W(t) = [)C(t), p(t)]+i[H(t), p(t)] + . We integrate numerically each matrix 
element of eq. (||) by the fourth-order Runge-Kutta method p| (thus we solve 
complex simultaneous differential equations). Wc use the following parameters: 
% = 1, 70 = 0.01,7i = 0.14,72 = 0.14, and A = 0.01. System size is taken to 
be 10 sites and the periodic boundary condition is imposed, that is, the system 
consists of two unit cells. All the numerical integrations are made with the time 
step of Runge-Kutta method to be 0.0001. The initial state of p{t) is prepared 
as the thermal equilibrium state of the inverse temperature (3. 

In Fig. ([l]) wc show a typical time series of the total current per unit cell 
3 = T, ln umtceii J n- Here we set f3 = 0.8 and T = 500. The first duration of 250 
time corresponds to the on-potential period in which the external force is zero, 
and next one for 250 time is the off potential period in which the external field 
cancels the asymmetric periodic potential. Oscillation in the current are clearly 
seen in Fig. ([!]): in the on-potential periods, fast oscillation appears but does 
not contribute to the time-averaged net current. Its typical time scale is the 
scale of h/AE, where AE is the typical eigenvalue difference in the asymmetric 
potential. On the other hand, much slower fluctuation is observed in the off- 
potential period. This is originated in the interaction Hamiltonian; its time scale 
is the order of I/70. In the off-potential periods, we recognize a synchronous 
appearance of the current, that is, every phase of the oscillations in each off- 
potential period has an almost same value. This synchronous oscillation of the 
current, which does not depend on initial conditions and has a positive value 
for the present parameters, mainly contributes to the net current when the 
long-time average is taken. 

In Figs. (|^) we show the time-averaged current (J), versus the inverse tem- 
perature j3 (a), and versus the period of the external field T (b), respectively. 
The current varies with (3 and T and takes negative or positive values according 
to the conditions. This result is in clear contrast with the classical thermal 
ratchets. In the latter, the direction of the current does not change with pa- 
rameters. In the present case the finite value of J results from the synchronous 
oscillation mentioned above. The direction and the magnitude of J depend on 
the phase of the oscillation at the instance the potential is turned off, which 
depends on the parameters. 

Contact with the heat bath is essential for the synchronous appearance of 
the current. In fact, such an synchronous oscillation disappear if we turn off the 
interaction between the ratchet system and the heat bath. Figure (||) shows the 
current as the function of time, as Fig. (^), but in the absence of interaction with 
the heat bath. In this situation, the ratchet system obeys pure quantum dynam- 
ics. The current still appears in the off-potential periods, but its magnitude and 
direction are random. The reason is the following: in the on-potential periods, 
the current oscillates according to the eigenvalue difference of the Hamiltonian. 
Since the period of external field and a oscillation period in the on-potential 
period is generally incommensurate with each other, the phase of the oscillation 
when the potential is turned off fluctuates almost randomly. 
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The mechanism of synchronous oscillation for the case of finite 70 is consid- 
ered as follows: After large enough repeated time n of the external field, the 
system becomes non-equilibrium steady state due to the dissipation. Then the 
density matrix p at each end of the period of the external field coincides with 
each other. It produces the synchronous oscillation. It is to be noted that the 
above stationary motion does not mean the stationary motion in each duration 
where the potential is constant. 

To summarize, we proposed a new concept, quantum ratchet, and its one 
possible realization. Finite net current is found in the non-equilibrium steady 
state, which originates from purely quantum mechanical effects. The behavior 
of the current is quite different from that in the classical one: in the quantum 
ratchet, the current flows both direction which is determined by the heat bath 
temperature and/or the period of the external field, though, in the classical 
situation, it flows to one-direction decided by asymmetry of potential and type 
of external field. The bidirectionality means that, in quantum ratchet devices, 
we can control the electric current by changing parameters ]23| . 

The present model has been considered in the weak coupling situation. Thus 
the relaxation time of the system becomes long. Because of this, we have studied 
the system only in the non-adiabatic time scale. If we study in the adiabatic 
and strong coupling region, different phenomena like crossover between classical 
and quantum regimes may be observed. The present results depend strongly on 
the detail of the model such as the potential shape, system size, and so on. This 
indicates that there is a possibility to produce the current more effectively by 
changing details of the model. For the study of the quantum thermodynamics 
of the rathcet system, especially properties as a quantum engine, or a quantum 
energy transducer, we have to study efficiency of energy [^4|. Detailed study is 
now in progress. 

The authors are grateful to S. Miyashita, T. Nagao, K.Saito for valuable 
discussions. This work is partly supported by a Research Fellowships of the 
Japan Society for the Promotion of Science for Young Scientists and also partly 
supported by Grants-in-Aid from the Ministry of Education, Science, Sports 
and Culture. 
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Figure 1: Typical time series of the total current per unit cell: initial state is 
chosen as a thermal equilibrium state. We take some parameters to be (3 = 0.8 
and T = 500. The system in first duration for 250 time is in "on-potential" 
state. In next one for 250 time it is in "off-potential" state. The total current 
J = ^2 in un i tce ii Jn hi the unit cell is calculated with the time interval 5. 



Figure 2: The average current (J) versus (a) the heat bath inverse temperature 
(3 and (b) the period of the external field T. In a classical overdamp limit, we 
expect only negative current to appear. 



Figure 3: Typical time series of the current under pure quantum dynamics, 
that is, in the absence of the interaction with the heat bath. Initial condition 
of the system density matrix is taken to be "thermal equilibrium" state with 
inverse temperature (3 — 0.8. All of the other conditions are same as Fig. (Q). 
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